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Abstract. We use the Real Space Renormalization Group (RSRG) method to study extreme 
value statistics for a variety of Brownian motions, free or constrained such as the Brownian bridge, 
excursion, meander and reflected bridge, recovering some standard results, and extending others. We 
apply the same method to compute the distribution of extrema of Bessel processes. We briefly show 
how the continuous time random walk (CTRW) corresponds to a non standard fixed point of the 
RSRG transformation. 



EVS from RSRG : Brownian Motion, Bessel Processes and Continuous Time Random Walks 



2 



1. Introduction 

There is recent a regain of interest for extreme value statistics (EVS) and its applications to physics 
[TJ [U [3J 21 [5] biology [7] and finance [HI W\ ■ Characterizing the statistical properties of the maximum 
I max (or the minimum X m [ n ) of a set of N random variables Xi, X2, • ■ • , Xjy has recently found many 
applications in these various areas of research, particularly in statistical physics. Indeed, being at the 
heart of optimization problems, EVS plays a major role in the study of disordered and glassy systems. 

The EVS of independent, or weakly correlated, and identically distributed random variables is 
now well understood, thanks to the identification, in the large N (thermodynamical) limit, of three 
distinct universality classes depending on the parent distribution of the X^s [10]. The difficult cases 
arise when the random variables are strongly correlated, which happens to be the case of interest in 
many problems encountered in statistical physics. This is the case for instance of the directed polymer 
in random media and its connections to Kardar-Parisi-Zhang equation [llj or the statistical physics 
of a single particle in a logarithmically correlated random potential [3J, closely related to Derrida's 
Random Energy Model [T2"] . 

The most important example of a set of strongly correlated random variables is of course Brownian 
motion (BM). And despite the fact that it has been widely studied [13l[T4l[T5], there has been a renewed 
interest in the study of functionals involving Brownian motions |16[ll7j. and in particular in the context 
of extreme statistics |T8j [19j [2Q1 EU [22] . Considering a Brownian motion u(x) on the interval x € [0, L], 
subject to u(0) = 0, it reaches its maximum u m at time x = x m . The basic question we ask is : 
what is the joint distribution PL(u m ,x m ) of these random variables? For instance, for unconstrained 
Brownian motion, the (marginal) distribution P[(i m ) is given by Pi(i m ) = ~Xm^ 2 {L — x m )~ 1 ^ 
or equivalently the cumulative distribution Proba(x m < x) — — sin -1 {^J 1 x / ' L) , which is the classical 
Levy's arcsine law [23] . These classical results for BM for PL(u m ,x m ) have been recently obtained 
using Feynman-Kac formula and further extended to a wider class of constrained Brownian motions, 
including Brownian excursions, meanders and reflected Brownian bridge 21 1. Interestingly, it was 
shown that this joint distribution Pi(w m , x m ) for a Brownian motion and Brownian bridge arises 
naturally in the study of the convex hull of planar Brownian motions [24] . Recently, such results for 
the distribution of the maximum Pl(m to ) have been obtained, using also a path-integral formalism, 
for the case of multidimensional processes, where one considers p non intersecting Brownian walkers 
(vicious walkers) [25[ [26] . 

On the other hand, it is interesting to extend these studies to other stochastic processes. A natural 
one is the Bessel process, which is the radius of d— dimensional Brownian motion. The study of the 
distribution of the maximum Pi(w m ) for a Bessel bridge has actually been widely studied in probability 
theory and statistics, in particular in the context of the generalization of the Kolmogorov-Smirnov test 
(corresponding to d = 1) 27J. The marginal distribution Pi(u m ) for such Bessel bridges was first 
computed by Gikhman and Kiefer 28, 29J and then generalized to any real d by Pitman and Yor [30] . 
while much less is known about Pi(i m ). In Ref. [28, 29} HQ], these results are obtained using rather 
complicated tools from probability theory, including the so called "agreement formulae" which, albeit 
rigorous, are not very physically intuitive (see however Ref. [21]). To our knowledge, the extreme 
statistics of the Bessel process has not been studied using the tools of statistical physics. 

In this paper we will use the Real Space Renormalization Group (RSRG) to study these problems 
of extreme value statistics. The RSRG was first devised for the study of quantum spin chains with 
disorder J3TJ [351 [33J [33]. It was then applied 35J to the study of the Sinai model [3S], i.e. the 
Arrhenius motion (diffusion) of a particle in a random potential u(x) which is itself a random walk, 
i.e. it becomes a one dimensional Brownian motion in the scaling limit. In both cases the idea is 
to decimate iteratively the smallest energy levels (for quantum systems) or barriers (for Arrhenius 
dynamics), defined as variables F = \u{x{) — u(xi+i)\ where (2^, £1+1) are pairs of successive of local 
min-max. At a given stage, indexed by the "decimation scale" called T, the landscape does not contain 
any barrier smaller than T. While it is only asymptotically exact for the above mentioned problems, 
it was later realized [37] that it is also a fully exact and quite powerful tool to study extreme value 
statistics, and that it leads to solvable equations when applied to a broad class of processes, which 
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includes, but is not restricted to, the Brownian motion. In fact it gives much more information than the 
frequently asked questions in extreme value statistics as it also allows one to compute the statistics of 
/i-extrema (as introduced in the mathematical literature by Neveu and Pitman 38J), i.e. the extrema 
separated by a barrier larger than h, and here h = Y. In Sinai's walk, these are "important extrema" 
of the energy landscape, which are separated by large enough barriers such that the process is not 
authorized to backtrack on more than T (see Ref. [39] for ft,-extrema in the context of Sinai's walk). In 
Ref. [37] the RSRG was used to obtain the extremal statistics of the toy model, i.e. Brownian motion 
plus a quadratic well. 

The aim of this paper is to explain and illustrate how the RSRG allows one to recover known results 
about extrema of the most standard constrained Brownian motions, i.e. the bridge, the meander, the 
reflected BM and the excursion. It is in a sense an overkill, as we will see the method is too powerful 
and we only use it in the limit T — > oo. However, in view of the versatility and importance of this 
method, we think that it is useful to establish some links with more classical approaches and results. 
Since we want to be mostly pedagogical, we restrict to symmetric BM, the case of a bias being a 
straightforward extension. Note that even for the Brownian motion we compute a rather complicated 
object, depicted in Fig. [21 at the same price as the standard, well known result for the distribution of 
the maximum. Next we study the Bessel process (i.e. the radius of a <i-dimensional Brownian motion) 
and obtain its extremal statistics. Finally we close by a short Section on the extensions to continuous 
time random walks (CTRW). Some of our results are standard well known results, derived here in a 
rather different fashion, others are extensions and likely to be not so well known. 

2. Extreme value statistics of some classical examples of constrained Brownian motion 

2.1. Short review of RSRG 

Let us briefly recall the main results of the RSRG method. Here we restrict to the Brownian motion. 
The general case, introduced in Ref. [37], is recalled later. The first statement is that the probability 
measure of the symmetric Brownian motion u(x) with x G [0, L] and u(0) = uq, u(L) — ul, can be 
written as the following sum: 

I r E r (F ll £ 1 )E r (F 2 ,£ 2 )S(L - (£ 1 + £ 2 ))5(u L - {u + F ± - F 2 )) 

oo 2fc+l 2fc+2 2fc+2 

+ Y / lrE r (F 1 ,£ 1 ) H P(F j ,£ j )E r (F 3h+ 2,ta h +. 2 )5(L- ]T £ t )S(u L - (u + £ (-l)^ 1 ^-)) . (1) 

k=l 2=2 i=l j = l 

where lr = Y 2 . In the (finite size) RSRG language !3j| this is called the renormalized landscape at scale 
r, with bond lengths £i — Xi~Xi-\ and barriers F 2n — u(x 2n ~i) — u(x 2n ) or F 2n+ \ — u(x 2n +\) — u(x 2n ) 
where x 2n are local minima and x 2n +i are local maxima, and it holds for a large class of random walk 
landscape with Markov property. It is illustrated in Fig. [T] The transformation which increases T is 
called the decimation. It iteratively reduces the list of extrema xq — < x\ < ... < x 2 k+ 2 < Xl = L 
by removing a pair of adjacent min-max (or max-min) from the list corresponding to the smallest F 
in the system, when Y crosses that value, in effect gluing three bonds into one. This results in a RG 
(i.e. evolution) equation for Pr(F, £) (bulk bonds) and Er(F, £) (edge bonds) which eventually reaches 
a fixed point form which, under some conditions, corresponds exactly to the Brownian motion. The 
edge bonds play a special role as they cannot be decimated, they can only grow by decimation of 
their nearest neighbor. The total measure ([T]) is a sum of measures for the events where there remain 
2k + 2 bonds in the system, i.e. k valleys of bulk bonds and two edge bonds, with k = 0, 1, ... At 
large Y, for L finite, there remains only the term k = which plays a special role. The product form 
reflects the Markovian nature of the landscape, which is preserved by decimation, hence is either an 
exact consequence of the choice of a Markovian initial landscape, or a consequence of the convergence 
to the fixed point landscape. The only constraint is the fixed total length, implemented by the delta 
functions. The factor lr = Y 2 is nothing but the average bond length which ensures the normalization 
of the total probability to unity. The fixed point form for the bond length probability, i.e. the one 
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Figure 1. Schematic representation of the finite size RSRG measure (TTJ for the renormalized 
landscape (on the x axis, the ij's are the bond length, and on the vertical y axis the Fj's are 
the barriers). It is a sum of terms with k = 1,2,.. valleys (i.e. sets of two consecutive bonds) not yet 
decimated at scale T, i.e. of barriers larger than Y (in vertical y direction). The edge bonds (first 
and last) cannot be decimated, only the bulk bonds (all the others) are. 



corresponding to the BM, takes in Laplace Pf(F,p) — e pi Pr{F,t), the simple form: 



Pr(F, P ) = 



VP 



sinh(Ty/p) 
E r (F,p) = r- 1 e- F V^oth(ry50 



e -(F-r)v?coth(rv?) 6)(_p_r) , 



(2) 
(3) 



As mentioned above Pr (F, t) can be related to a sum over all BM which are constrained to climb from 
to F in a time i while remaining in the interval ]0, F[ and with no descent of more than T. The 
RSRG equations which yield the fixed point forms [3J are recalled in Appendix A. 



2.2. Free Brownian and Brownian Bridge 

In all cases one considers a Brownian motion u x on the interval [0, L], with (du%.) = 2dx, with initial 
condition uo = 0. Note that this corresponds to a diffusion coefficient D = 1, in contrast to the 
value D — 1/2 used in Ref. [H] to which we will refer in the following. While the free Brownian is 
unconstrained at x = L (see Fig. [2|, the Brownian Bridge is defined by the constraint uj, — (see Fig. 
fj]a)). It is natural to consider the joint probability density PL(u m ,x m , Ul) that the free Brownian 
ends at u l and has maximum value u m at position x m . From the above considerations it is obtained 
from RSRG by considering the last block (see Fig. [3]), and is given by: 

P i (ii m ,x m ,H I )= lim T 2 E T {F x ,x m )E T {F 2 ,L-x m ) , (4) 

with F\ = u m and F2 = u rn — uj^. Using the result ^ for the edge bond, one finds, in the large T 
limit: 



lim TE T (F,p) = e- F ^ , 



(5) 
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Figure 2. Unconstrained Brownian motion u x starting at uo = and ending at ill which is free. It 
reaches its maximum u m at x m and has a minimum ui, reached at xi, on the segment [0,x m ] and 
a minimum U2 , reached at X2 on the segment \x m , L] . 

which gives after inverse Laplace transformation : 




X m L X 

Figure 3. Schematic representation of the renormalized landscape associated to the one in Fig. (f2j) 
after the last RSRG step, i.e. the last decimation. The measure of this last block gives the joint 
distribution P(x m , u m , u^) (see Eq. ©). 

Thus using Eq. (f4|) together with the expression for the edge bond in Eq. (J6j) one has 

P L (u m ,x m ,u L ) = ^ '—e 4 — *(i-»m) 0(u m - u L )9{u m ) . (7) 

4?r x l( 2 (L-x m ) 3 / 2 

Integration over x m and u m yields the propagator W(ul, L) 

W(u L ,L)= du m dx m P L (u m ,x m ,u L ) = —==e~^ . (8) 

./max (0, tti,) JO V47rL 
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We now obtain the result for the free Brownian motion by integrating the formula ([7]) over the 
final position ul- 

f Um Urn "™ 

Pl (u m ,x m )= / du L P L {u m ,x m ,u L ) = 9(u m ) m 3 / ^ e ~ I ^" r ■ ( 9 ) 

J — oo 2tT\/L — X m Xm 

This yields the well known results for the marginal distribution of the position of the maximum x m 
for the free Brownian motion 

VT c {x m ) = \v b ™ fe) , f hcc (x) = 1 , (10) 

\ Jj / 7TWX{1 — X) 

which is the well known Levy's "arcsine law" [23j , as well as the distribution of the maximum u m 

x 2 

p l r c M = ^p fr - (^) , p ft ™(x) = e(xf-^- , (ii) 

where 9{x) is the Heaviside step function, 0(x) = if x < and 9(x) = 1 if x > 0. 

From Eq. ([7]) we also immediately get the result for the Brownian Bridge, where ul = (see Fig. 
Ha)): 

P? idsc (u m ,x m ) = P L (u m ,x m ,u L = 0) = JT_ e __^_ 



(12) 



This yields the marginal distributions of a; m on the one hand and of u m on the other hand 



V h L Mgc (x m ) = \ , (13) 



pbr 



idgc (« m ) = j^P hlldse (^) , P brids °(^) - 20( 2 ;)^e- a:2 . (14) 

We thus see that the RSRG method allows to recover standard results for the relatively easier cases 
of free Brownian motion and the Brownian bridge. Let us now investigate more complicated instances 
of other classical constrained Brownian motions. 

Notational remark : in the following we will use invariably the notation Pl{- ■ ■), independently 
of the number of its arguments, to denote some joint distribution. Similarly, we will also use the same 
notation to denote its Laplace transform with respect of one or several of its arguments. The choice 
of the notations for the different variables should be clear enough to avoid any confusion. 

2.3. Reflected Brownian motion, Brownian excursion and meander 

In this section, we study the cases of the reflected Brownian motion, excursion and meander, see Fig. [4] 
b), c) and d) respectively. The joint distribution PL(u m , x m ) of the maximum u m and its position x m 
for these three cases was recently computed using both path integral techniques as well as the so called 
rigorous "agreement formulae" in Ref. [21]. Using RSRG, all these cases can be obtained from the 
probability Pt(ui, x±, u m , x m , u 2 , x 2 , ul) that an unconstrained Brownian u x starting at uq = ends 
at ul, reaches its maximum u m at x m and that its minimum U\ on the segment [0,x m ] is reached at 
X\, and that its minimum u 2 on the segment [x m , L] is reached at X2 (see Fig. [5]). It can be obtained 
from the edge bond probabilities (including information on the minima) as explained in |Appcndix C| 
as: 

P L (ui,x 1 ,u m ,x m ,u 2 ,x 2 ,u L ) = E(u m ,x m , -Ui,Xi)E(u m -u L ,L- x m , -u 2 + U L ,L- x 2 ) , (15) 

lim TE T (Um,Xm, u,x) . 
We show in |Appcndix C that in Laplace variables: 



poo poo 

E(u m ,p,u,q) = / dx m / dxe~ pXm ~ qx E(u m ,x m ,u, x) , 
Jo Jo 

BV \ sinh(yp + qu rn ) y/p 

E{u m ,P,u,q) = . /= =} — ; ■ rr n , — ; \T • (16) 

smh( % /p + q{u + u m )) smh(^/p(u + u m )) 
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Performing Laplace inversions (see Eq. (jC.121 |C.14[) ) one finds: 

E(u m ,x m ,u,x) = J2 4 7 r 3 nm 2 (-l)"+- 7 ( ^5 ) e "^ 7( " W(Xm ' Z)) - ( 17 ) 
Notice the scaling forms, directly read from Eqs. (|16[ [T7|) 

E(u m ,x m ,u, x) = —E ( — , j , (18) 

1 ~ / , U 

E{u mi p, u : q) = E pu m , , qu ri 



U 7n \ U r 
c x 



m I ' 



poo poo 

with E(p, u, q) = / dx m / dxe~^ Xm ~^E(x mi u, x) . 
./o Jo 

From Eqs. (fT5l fT7l) one gets the probability density of the event depicted in Fig. [2J 

oo 

P L (u 1 ,x 1 ,u m ,x m ,u 2 ,x 2 ,u L ) = 6>(u L - u 2 )0(u m - u £ ) l&n & nin 2 m\m 2 2 {-l) ni+mi+n 

fii ,"^1^2 ,m2— 

x gl^^g^^ (ng(£-» 3 )+m3(s 3 -* m )) (1Q) 

(U m - Ux^iUjn - U 2 ) 5 

From this object (|19p we now obtain the probability density PL(u m , x m , u, x, ul) that an unconstrained 
Brownian u x starting at uq — ends at u £ , reaches its maximum u m at x m and its minimum u at a;. 
It is given by 

P L (u m ,x m ,u 7 x,u L ) = 0(x m - x) du 2 dx 2 P L (u, X, u m , x m , u 2 , x 2 , U L ) (20) 

J u J x m 

+ 0(x~x m ) / dux / dxiP L {ui,xi,u m ,x m ,u,x,u L ) . (21) 
After some algebra, left in | Appendix D| one obtains 



P L (u m ,x m ,u, Xl u L ) = 0(x m - x) 



(u m — u) 5 



ri2— 



■ TT(u L -u) n ^ ^2^2 

{u m - u) 2 



x V 27r(-l)" 2+1 n 2 ; ^-n^, ^——^^-^) 

^— ' ill — 



+ 9{x-x m ) £ ^n 2 ml (- 1)"^ S 'f " 2 } e ~ ("I (x-x)+r^(x- 8m )) 

712, m, 2 

x V 2 7 r(-l)" 1+1 n 1 9 iy e -(^F a '" . (22) 

ni=0 x ' 

In particular, one can check explicitly the following symmetry 

P L (u m ,x m ,u,x,u L = 0) = P L (-u,x,-u m ,x m ,u L = 0) . (23) 

From the structure of Eq. (|15p one obtains the useful expression for -Pj,(lt TO , x m , ul, a) which is 
the joint probability of the maximum u m and of its position x m for a Brownian motion starting at 
uq = and ending at ul and with a minimum value u larger than —a. It reads : 

du I dxPL(u 

-a JO 

ul+cl pL—x m pa rXm 

du / dxE(um — ul, L — x m ,u,x) / dv! I (ix'£(u m ,x m ,w',i') . (24) 
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Figure 4. Various constrained Brownian motions considered in this paper : a) Brownian bridge 
which is a Brownian motion constrained to start and end in 0, i.e. which no = «l = 0. b) Reflected 
Brownian bridge, u x = where u® is the standard Brownian motion, and constrained to start and 
end in 0. c) Brownian excursion, i.e. Brownian motion constrained to start and end in 0, remaining 
positive in between, d) Brownian meander, Brownian motion starting in 0, ending anywhere in 
Hi > and remaining positive in between. 



2.3.1. Reflected Brownian bridge The reflected Brownian is defined hy u x = where u^. is the 
usual Brownian, and initial condition uq — 0. We focus on the reflected Brownian bridge, for which 
uj, = 0, see Fig. S]b). Therefore one has 



P™ ' " se (u m ,x m ) — V4nL( / du I dxP L (u m ,x m ,u,x,u L = 0) 



o 



dUm I dx m Pl J (u rni X m , W m ,X m ,Zi£ — 0)j , 

o Jo 



(25) 



(26) 



where the prefactor v AttL comes from the fact we are considering bridges (similarly to Eq. (fT2|) above). 
Using the above symmetry ([23]) one sees that these two terms are actually identical, yielding 



P™ ' n ee {u m ,x m ) — 2V47r£ / du J dxP L (u m , x m , u, x, u L = 0) . (27) 



Therefore one has P, 



ref .bridge , 



u m ,x m ) — 2\/47rLr\(u m , x m , uj, — 0,u m ), and the above identity in 



Eq. ((241 y ields 



P^ e ' n ge (M m , x m ) = 2V AttL I du 



L-a 



dxE(u m , L — x m ,u,x) / du dxE(u m ,x m ,u,x) . 
Jo Jo 

(28) 

From the scaling forms in Eq. (fT5|) one obtains that 



^ref . bridge 



(u m ,x m ) = 2V4^LP C )P C 



L — u Tl 



(29) 
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where the Laplace transform of P c (z) is given by 

o *s-P M . [ d uE {p , „ . 0) . £ dUj -^My . _L_ . (30) 
By inverting the Laplace transform, one obtains 

oo ^ 

P c {z) = 7T ^(-l) n (n + -) e - 7r2 ("+ 1 / 2 ) 2z . (31) 
And finally, one has 

D ref. bridge/ % , 5/2 r 1/2 1 / ■wn+mi , ^ \ I , ^ \ -ff 2 (m+l/2) 2 ijp- -jr 2 (n+l/2) 2 

Um m,n=0 

(32) 

which gives back the formula obtained by Majumdar et al. in Ref. [3T] (see their formula (28)). In 
particular, by integrating over u m , one obtains the marginal distribution of x m as 

(z m ) = ^ ref - bridse (^) , (33) 

prof. bridge = 2 y ( _ ir +« (2m+l)(2n+l) _ 

^tlo [(2n + lfx + (2m + 1)2(1 - *)] 3/2 

2.3.2. Brownian excursion Here one considers the case of Brownian excursions, constrained to start 
and end in 0, remaining positive in between, see Fig. |4]c). In that case, one would naively impose 
that the minimum is zero exactly and ul — 0. However, one immediately sees on Eq. (1201) that 
PL{ u m,x m ,u = 0, x — 0,it£ = 0) = 0. This is expected since the Brownian motion has an 
infinite density of zero crossings. Therefore, to compute the joint probability P^{u m ,x m ) from 
PL(u m ,x m ,u,x,UL = 0) in Eq. ([2H)l . one uses a limiting procedure (see Ref. [T5J, [2T] for a similar 
procedure used within Feynman-Kac formula): 

pcx, \ r fe du Jo dxP L { Um , X m ,U, X, U L = 0) 

"£ \u m ,x m ) — am — ^ =j — g ^ . (35) 

£ ~ +0 f™du m J dx rn J ( du /„ dxP L (u m , X m , U, X , U L = 0) 

Following the same analysis as above, and using again the identity in Eq. (|24p one obtains that 

du I dxP L (u m ,x m ,u,x,u L — 0) (36) 



•rjrcf . bridge ' — 
7 L 



JO 

pL—x m rXm 

e 2 dxE(u m , L - x m , 0, x) / dxE(u m , x m , 0, x) 4- C(e 3 ) 

Jo Jo 



X 7f 



Using the above scaling forms (fT5)) one obtains that 

du dxP L (u m ,x m ,u,x,u L = 0) = £ 2 -^-P s ( j P 4 

Je JO M m \ X mJ 

where the Laplace transform of P s (z) is given by 

dze-^P s {z) = u = 0, g = 0) = — • (38) 
By inverting the Laplace transform, one obtains 

oo 

P s (z) = 2tt 2 {-l) m+1 m 2 e- m2 * 2z . (39) 

m— 

Finally, after normalization, one obtains 

3^9/2^3/2 - 



PI x (u m ,x m ) = °" £ ]T (-l) m+n m 2 n 2 e- m * ^e"" " ^3T , (40) 
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which gives back the formula obtained by Majumdar et al. in Ref. [ST] (see their formula (8)). In 
particular, by integrating over u m , one obtains the marginal distribution of x m as 

VT( Xm ) = ^(^) , (41) 

00 2 2 

V cx (x) = 3 V (-1)'"+" — — . (42) 

„,m=o [n 2 x + m 2 (l~ x)] 1 

2.3.3. Meander Here we consider Brownian meanders, which are similar to excursions except that in 
that case the endpoint ul is free, see Fig. |4]d). Thus one would impose the minimum to be and 
integrate over u^. Following the same limiting procedure as above one has 

Pl™(u m , Xm ) = Um - g dUL £ d f L ^ X ~ UL) • (43) 

J" m du L J °° du m f Q dx m J e du j Q dxP L (u m ,x m ,u,x,u L ) 

Performing the same analysis as before, and using again the identity in Eq. (|24p one obtains that 
du L du dxP L (u m ,x m ,u,x,u L ) 



o 

ru m rL—x m r^m 

= e du L / dxE(u m - u L , L - x m , 0, x) / dxE(u m , x m , 0, x) + C(e 3 ) . (44) 
Jo Jo Jo 

Therefore, after straightforward algebra, one has here 

1 „ ( Um \ _ / L u r , 



du L du dxPL(u m ,x m ,u,x,UL) = £ — -Ps I — I Pt I j ) ' ( 45 ) 

Je Jo u m \ x mJ \ x rn / 

where the Laplace transform of Pt{z) is given by 

r 1 r du tanh(^) 

Inverting the Laplace transform yields finally 

4_5/2rl/2 - 2 2x m 2 2t-. m 

P?-(«m,*m) = E ((- 1 ) m+ " - (-l)")^ 2 e-" * " T , (47) 

,ft m,n— 1 



which gives back the formula obtained by Majumdar et al. in Ref. [2JJ (see their formula (19)). In 
particular, by integrating over u m , one obtains the marginal distribution of x m as 

x I. 

~L 

,2 



Vr a M = V™ a {^) , (48) 



n 

[n 2 x + m 2 (l — x) 

,2 



V™*{x)= E ((-!)-+" -(-!)") — - (49) 



= 2 V (-1)" +1 - — , (50) 

m=o£=i Kx+ (2m + 1)2(1 -x)] 3/2 

where in the last equation we have used that only the terms where m is odd contribute, i.e 
(— l) m +" — (—1)™ = when m is even. 

We have thus shown that the RSRG allows one to recover in a rather simple way standard results 
for various constrained Brownian motions. We now extend this method to other stochastic processes 
and consider the case of Bessel Processes. 
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3. Bessel Processes and more 



3.1. RSRG for more general process 

The RSRG was extended to more general processes in Ref. [37], and we refer to that paper for all 
details. The generalization of the finite size measure (JT|) now reads: 

E^(0,u ;ui,Xi)E^(ui,Xi,L,u l ) + (51) 

OO k 

-Er (°i u o; ui, xi) JJ -Bjt" {u2j-i, X2j-i; u 2j , x 2 j )S r : (u 2 j , x 2 j ; u 2 j+i, x 2j+1 )E^ (u 2k+1 , x 2k +i , L, u L ) , 
k=i j=i 

which has still a block structure, but with slightly more general blocks. The index "— " refers to an 
ascending bond, i.e. where u x is on average an increasing function, while "+" refers to descending 
bonds. This distinction was not needed and suppressed in |T]) for the symmetric BM, but it is useful in 
general. As compared to Ref. [37], we have reversed the orientation ±. This amounts to the reflection 
u{x) — > — u{x) w.r.t. the notations there, or, in other words to consider here u(0~) = u(L + ) — — oo. 
This implies that the k = term (the first one) in the large T limit yields the distribution of the 
maximum (rather than the minimum there). RSRG equations for these blocks were derived there, 
which we will not reproduce here, and only a small class of solutions were identified (presumably many 
more remain to be discovered). 

As an example, the following class of real valued Langevin process u x was found to provide 
solutions of the RSRG equations, and was studied in Ref. [37] : 

du x = - W'{u x ) + dB x , (52) 



with dB 2 — 2dx a Brownian motion. It represents processes which undergo diffusion in a given 
potential W(u) with u s] — oo, +oo[. The explicit form of the blocks, for any T, was obtained in Ref. 
[37], and we need here only the form of the edge blocks, given in Laplace w.r.t. x\ — xq as: 

E±(u , Ul ,p) = e-i( w ^- w ^E±(u , Ul ,p) , (53) 
Ey (uq, u,p) = E^(u, uq,p) = exp (^J dvd\ \nK(v, v — T,p) 



where 

K(u 7 v,p) = — ^ ($i( M ,p)$ 2 (w,p) - $ 2 (u,p)$i(v,p)) , (54) 
where $i(u,p) are the two independent solutions of the associated Schrodinger problem: 

dl-(p + V(u))Wu,p) = , V{u) = ^W'(u) 2 - )^W"(u) , (55) 

and w(p) their Wronskian, w(p) — < &i(u, p)d u & 2 (u, p) — <& 2 (u,p)d u <S>i{u 1 p). Bulk bonds have a similar 
expression in terms of the kernel K(u, v,p). Note that in the limit Y — > oo the quantity: 

di mK(v,v - T,p) = — — — — - — — — — — -, (56) 

®i(v,p)<!> 2 (v-r,p)-<f> 2 (v,p)<f>i(v-T,p) 

automatically becomes equal to 3>'(v,p)/$(v,p) where <5>(v,p) is the linear combination of $i and 
$2 orthogonal to the one which blows up at v — > — oo. This is what is expected for a process 
u x g] — oo, +oo [. 



3. 2. Extreme value statistics of the Bessel Process 

We now call u x the Bessel process, i.e. the radius of the c?-dimensional Brownian motion u x = 
\f}2i=i{B x ) 2 where the B l x are d independent BM with (dB l x ) 2 = 2dx, i = 1, . . . ,d. As is well 
known it satisfies the one dimensional Langevin equation du x = dB x + (d — Y)/u x , u x > 0, with 
dB^ = 2dx. Hence it is of the type ([52"]) with a potential W(u) = — r In \u\ with r = d— 1. This allows 
a generalization for real values of d. 
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Figure 5. Bessel process in dimension d = 2, i.e. v = 0. a) : Brownian motion in dimension d = 2. 
b) The Bessel process with index v = corresponds to the radius of the 2-dimensional Brownian 
motion. 



Consider now a Bessel process (BP) starting at u = uq at x = 0, and ending at ml at x = L. Let 
us call Pi(M m , £ m , Ui|uo) the joint probability distribution that the BP u x starting at uq has minimum 
u m at position x m and ends up at ul ■ It is given by the last RSRG block in (|51[) : 



P L (u m ,x m ,u L \u ) = lim E r (u ,u m ,x m )E^(u m ,u L ,L - x m ) . (57) 
r^oo 

It satisfies the normalization condition: 

rL ru M 

dx m du m du L P L (u m ,x m ,u L \u Q ) = 1 , (58) 

J — oo J it m 

which we will check below explicitly on our solution. In the case of the BP, there is a subtlety with 
applying the recipe given in Ref. [37] and the solution (|5"5|) . Indeed the process remains on the interval 
u x E [0, +oo[. Hence u — has a special role, related to the returns to the origin. 

Let us first examine the naive solutions of l[55p for the potential V(u) = jW'(u) 2 — jW"(u) = 
(d — l)(d — 3)/(4m 2 ) for the BP. The two independent solutions read: 

$i(u,p) = y / uK 1 _d(u^/p) , $ 2 (u,p) = Vul 1 _d (uy/p) , (59) 

where I p {x) and K p (x) are modified Bessel functions of the first and second kind respectively. Note that 
for d = 1 one recovers the result for the symmetric BM, in which case W = and $i(u,p) ~ e~ u ^, 
< I > 2 {u,p) ~ sinh(u^/p), up to unimportant p-dependent normalization which is cancelled by the 
Wronskian. For d < 2 the above solutions behave as u — > and for fixed p > 0: 

$i(«,p) = 2- ,i/2 j)T 2 iiTr(l-rf/2) , $ 2 («,j)) = 2- 1+li/2 ])T ;! u ¥/r(2-d/2), (60) 

which, upon a rescaling by p( d_1 )/ 2 gives also the p = limit. In the above expression, T(x) is the 
Gamma function. In that limit p = these two solutions correspond to the general expression 37J: 

$i(«,p = 0) =e-s w (") , $ 2 (u,P= 0) = e-i ff W / du'e^"') . (61) 



The process being only defined for u x > the question arises of the correct solution which is 
selected in the limit T — > +oo. First of all the above formulae in Eq. ([53]) must be replaced by: 



E±(u , Ul ,p) = e-i( w ^- w( - u °»E±(u , Ul ,p) , (62) 
i?P (uo, u,p) — E^ (u, Ua,p) — exp / lnK(v, max(0, v — T),p) I , (63) 



and the question of the value for K (u, 0,p) arises. One way to regularize the problem would be to add 
near u > a very steep barrier (e.g. W(u) — e u °^ u ^) and consider K(u,v,p) in the limit \v\ 3> u 
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both going to zero. A short-cut to treat that barrier problem is to notice that this is equivalent to 
forget the barrier and ask for reflecting boundary conditions at v — for the process u x . We must also 
ask that as u — ► the current vanish so that: 

J s = (d u - hw'iu^J K(u,O lP ) = (d u - K(u 7 :P ) = . (64) 

It is then easy to find that the proper solution is K(u,0,p) ~ cosh(u^/p) for the symmetric BM (i.e. 
when W(u) = i.e. the d = 1 BP becomes the reflected BM see below), while for the BP we find: 

iir(«,O,p) = A0Bin(d|)$i(«) + *2(u)J = <M-i+i {uVp) ■ (65) 

where A is some constant. We may understand this from a limit case of the regularization by a barrier 
at scale uq writing: 

K(u,v,p) ~ $?°(u)^°(u) - ®f(u)$X°{v) , (66) 

where $"° 2 (u) are respectively decaying and exploding at u = +oo, and showing that irrespective of 
the form of the barrier it reflects hence: 



lim 



Using this result we now find that the large T limit is given by: 

2 „• fird 



-sin (67) 



lim Er(u ,u m , P ) = e -i^)-^» » ffiM + 



r^oo 



~ sin(df )$i(u m ,p) + $ 2 («m,p) 



d-2 2 



3in(^)JiT 1 _|(Mo v ^) + I 1 _*(u 0y /p) 
u 7 f sin( 2 ^)l£' 1 _|(u mv /p) +/ 1 .|(« mv ^) 



"0 / I- 1+ $(u m y/p) 

which is found to hold for any d. One finds similarly: 



(68) 
(69) 
(70) 



lim E^(u m ,u L iP) = — ^ — . (71) 



One can check the normalization (|58|) . i.e. in Laplace, i.e. one writes the product: 

r hm £t K<w)^ (wl.p) = -^)75 /_ 1+f(MmV p )2 ' ( 72 ) 

Integration over it£ and described in (|55|) yields 1/p, as required. 

For d = 1, which is just the reflected Brownian (see Fig. [4]b) one finds: 

r t?-( \ cosh ( M "\/P) r n cosh(u L% /p) 

hm £ r u ,Wm,p = T7 ft i }™. E^(u m ,u L ,p) = — — . (73) 

r^oo cosh(u mv /p) r^oo cosh(w mv /p) 

Similarly for d — 3, one finds: 

,. jp-f s u rn sinh(M 0v ^) . m £ sinh(u L ^/p) 

hm E r (u ,u m ,P) = r-TT pr , J™ E^(u m ,u L ,p) = — — . (74) 

r—oo u smh(M mv /p) r^oo u m smh(u mv /p) 

The Inverse Laplace transform now gives, for any d: 

lmX o E r (u ,u m ,x) = ^—j IUihn) e " m > W 
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with I v {z) = e~^™ v J v {iz), v = — 1 + d/2, and < j„ t i < j v .2 < ... is the sequence of positive zeroes 
of J„, i.e. Ju(jv,n) = 0, where J u {x) is a Bessel function of the first kind. It can be further simplified 
using standard relations between Bessel functions 

V JP-I \ f u m\~ST^ 2 ■ Jv( jv,n^) , , 

lim E r (u ,u m ,x) = — Z^~T^^ n ~i T- — \ e > ( 76 ) 



and, similarly: 



lim Ei(u m ,u L ,x) = ( — ) V" -|~jv,n T^ U ? Um \ e "™ ■ 1 ' ' 1 



Combining these two expressions (|76[ I77|) yields our final result for the joint distribution of the 
maximum, its position and the endpoint of the Bessel Process: 

pi I \ 4 • J v\Jv,n Um )-Jy\Jv.m Um ) t-, •— i 

mu m ,Xm,U L \u Q ) - ^_ 2)/2 ^ Jv+i y Wfn ) Jl)+1 ^ iTO ) 6 m 6 " < 

(78) 

with z/ = -1 + d/2. 

We now focus on Bessel bridges for which uq = = 0. Here again one uses a limiting procedure 
as before and writes 

pBessel bridge^ ^ = _ P L (u m , X m ,U L = e\u Q = e) (?g) 

/ e °° rfum Jo dx rn P L (u m , x m , u L = e\u Q = e) 
Using the small x behavior J v (x) ~ x v j(2 v Y(y + 1)) one obtains 



jBcsscl bridge/ \ \ * Jiak J^.m ^ 

(Wm,a;m) " r(i + J„ +1 (j„, n )J„ +1 (j„, m ) e 



Pr Se DTiaBe (um,x m ) = -— ^— „ > t "'Yt V , e ^e . (80) 



For d = 1, one has ^ = —1/2 and j-x/2,n = (n + 5) 71 ": one checks that this formula gives back the 
result for the reflected Brownian motion in Eq. (|32[) . For d = 3 one has v = 1/2 and j—x/2,n = n7T ' 
one checks that this formula gives back the result for the excursion in Eq. (|4T)]) . 

On the other hand, if one integrates over x m one recovers the distribution of u rn as obtained by 
Gikhman [28] and Kiefer [29] for integer values of the dimension d and later generalized by Pitman 
and Yor in Ref. [30] for non integer d (see |Appendix E[ ) . It reads 

^^<*>-z(iF^gdb"»H?)) ■ (81) 

Now if one integrates instead over u m one obtains the marginal distribution of x m as 



V 



Bessel bridge / 



(Z m ) = i^Bcsscl bridge /^j (82) 

L V L J 

^Bessel bridge^) = ^ + „) £ . j^j^m _ _ _ ^ 

n,m Jv+x(jv,n)Jv+x{jv,m) \_3v,n x ~T" jv,m(^ ~ ■")] 

Many more quantities could be computed and it remains also to extend these considerations to 
finite T and bulk bonds. We leave this to future studies. 



4. Continuous time random walks: a new RSRG fixed point 

We now come to the study of the so called continuous time random walks (CTRW), which was first 
introduced by Montroll and Weiss in Ref. [42] . Within this model of diffusion, the walker performs a 
usual random walk but has to wait for a certain "trapping time" r before each jump. The trapping 
times between each jump are independent and identically distributed random variables with a common 
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density function i/}(t) which has a power law tail, ip{ T ) ^ T~ l ~ a with < a < 1, while the jumps 
themselves are distributed according to a narrow distribution. This type of model was suggested by 
Scher and Montroll [43] to model non-Gaussian transport of electrons in disordered systems and since it 
has been widely used to describe phenomenologically anomalous dynamics in various complex systems 
[40 \ I41j . Indeed, for a < 1, the mean trapping time between two successive jumps is infinite and 
hence CTRW is characterized by a subdiffusive behavior, with a dynamical exponent z — 2/ a > 1, 
and non-Gaussian statistics. While CTRW has been widely studied [101 HI], it seems that the extreme 
statistics of these processes have not been studied in detail. It is the purpose of this section to study 
them using RSRG. 

We now trade the variable r for x. The CTRW described above generates a landscape u(x). 
Performing the first stages of the RSRG method as described in Ref. [37] , one finds that the bonds in 
the renormalized landscape at scale T acquire a broad distribution, with a typical length T 2 / Q , while 
the jumps in u remain of order T. One can now search for a fixed point of the RSRG which has 
precisely this scaling with T. This is done in Appendix B| and we find that these new fixed points can 
be simply obtained (up to some subtelties explained below) from the Brownian case by substituting p 
to p a (a = 1 corresponding to Brownian motion) in the RSRG-blocks of the BM. 

Therefore it is natural to consider the joint probability density QL(u m , x m , «l) that the free CTRW 
just arrives at ul at time L and has maximum value u m at position x m . From similar considerations 
as in Eq. ((4]), and the discussion in | Appendix B[ it is given by 

Q L {u m ,x m ,u L ) = lim ry a Er(F u x m )E r {F 2 ,L-x m ) , (84) 

r— >oo 

with F± — u rn and F2 — u m — UL- Using the result ([3]) for the edge blocks together with the substitution 
p to p a (see also Eq. IB. 61 in |Appcndix"Fj| one obtains in the large T limit: 

lim T^ a E r (F,p) = e- Fpa/2 , (85) 
r^oo 

which gives after inverse Laplace transformation : 

lim T^ a Er(F,l) = r 2 /"4 /2 (ir 2 / a ) , (86) 
r^oo 

where C v {x) is a one-sided Levy stable pdf of index v — a/2, whose Laplace transform is 
Jo°° e ~ px £v{x) = e~ pU . The function C v {x) can be represented in terms of Fox H functions [44l [45l 146] 
or, alternatively, can be written as a power series [47j 

C v (x) = - Y(-i)"+i x -™-i T ( n " + sin {rmv) ,x>0, (87) 

n=l 

from which one obtains the asymptotic behavior for large argument x ^> 1 as C v {x) ~ vT(l — 
v)~ 1 x~( 1+u \ where we have used r(^)r(l — v) = it/ sin (vtt). For small argument, it behaves as [44] 

C u (x) - x- a e- KX ~ T , (88) 

where r = vj{\ - u), k = (1 - v)v v /^- y \ a = (2 - i/)/(2(l - v)). 
For v — 1/2, one has 

W) = v^ e - 1/(4 ^ . ( 89 ) 

and thus the formula ([56]) yields back the result for Brownian motion in Eq. ©. For other values of v, 
a representation of C u (x) in terms of other simpler special functions is usually obtained by examining 
the series expansion in Eq. (|87p and using properties of gamma- functions. For any rational value of 
is, it was shown in Ref. |43] that C v (x) can be expressed in terms of a finite sum of hypergeometric 
functions (for instance, C\u{x) is given in Ref. [46]). Here we simply quote another interesting simple 
case corresponding to v — 1/3, i.e. a = 2/3 for which one has 

A/a(x) = 7^)T7J Ai [(3*r 1/3 ] , (90) 
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where Ai(x) is the Airy function. 

Now the probability that the free CTRW ends in «i at time L and has maximum value u m at 
position x m is obtained from QL(u mi x m ,UL) as a convolution : 

PL{u m ,X m ,UL) = / dL'QL'{u m ,Xm,UL)^{L - ll) , (91) 



where ^(l) — tl)(x')dx' is the probability of no jump in the time interval [0, 1] . Its Laplace transform 
^(p) behaves, for small p, as ^(p) — Ap a ~ x with A = 1, see | Appendix B| And therefore one obtains 
the double Laplace transform of Pi,(u m , x m , Ul) with respect to x m and L as 

/>L />oo 

dx m e-^ / dLe- pL P L (um,a; m ,«L) = e -(P+^«™ e -^("™— V^ 1 , (92) 



where v = a/2. From Eq. ([92| . one checks that one recovers the correct expression for the propagator 
W{ul, L) 7 namely the probability that the free CTRW ends in ul at time L. It is given by 

/>oo />L 

W(u L ,L)= / du m / dx m P L (u m ,x m ,u L ) . (93) 



./ max(«t ,0) J 

One obtains its Laplace transform with respect to L as 



W(ul,p) - /°° dLe~ pL W(u L ,L) = lp^ e ~P^ , (94) 
Jo z 

which yields the correct expression for the Laplace transform (with respect to time) of the propagator 

[4Tj . One can invert the Laplace transform to obtain 

W{u L ,L) = —J——£„ (—^—) ,v = a/2. (95) 



For a = 1, it yields back the expression in Eq. ([5]) and for a = 2/3 Eq. (|95p takes the simple form 

^K^) = (JgL) , (96) 



2iV3 V(3i) 1/3 , 

with Ai(x) - (3 2 / 3 F(2/3))- 1 - x(3 1 / 3 r(l/3))- 1 for z -» and Ai(x) - (2a; 1 /2 7r i/2)-i e -2x 2 / 3 /3 for 
a; — > oo. 

We consider the case of the free CTRW where one integrates over the final position ul between 
— oo and u m . In that case one immediately gets, from Eq. (|92|) the expression for the joint distribution 

P^ cc (^H m ^ Xrri) &S 

Pt c (u m ,x m )=8(u m )8(L-x m ) 1 1 -^fe) • (97) 

In particular, from Eq. (|97p one obtains the marginal distribution of x as 



V f rM = y^ frcc (^) , V^(x) = 1 , (98) 

L V L J ir x v m "(1 - x m )^ 

where we have used Jq 00 dyy~" ' C v {y) = 1/T[1 + i/]. Similarly, one has the marginal distribution of the 
maximum u m as 

Pl°>™) = ^;P froc (^) , P frcc (^) - e{x)±x-( 1+1 ^L u (x- 1 ^) . (99) 

Of course, following the same lines as done previously for Brownian motion, one could also consider 
various constrained CTRW: this is left for future investigations. 
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5. Conclusion 

In conclusion, we have shown that the RSRG is a quite powerful method to compute the extreme 
statistics of various physically relevant stochastic processes. By exploiting the solution of the RSRG 
equations found in Ref. |35H37j . we have studied the extreme statistics of the one-dimensional Brownian 
motion (BM) as well as the Bessel process (BP), i.e. the radius of the d-dimensional Brownian motion. 
For the BM, we have shown that it allows to recover, in a rather different way, standard results for 
the Brownian motion and its variants, including Brownian bridge, excursion, meander as well as for 
the reflected Brownian motion. For the BP, we have recovered in a simpler way the results of Pitman 
and Yor for the distribution of the maximum, and obtained also the distribution of its position. 

We have then extended this method to study the extreme statistics of the Continuous Time 
Random Walk (CTRW), which we have shown to correspond to a new fixed point of the RSRG 
transformation. Although we have restricted our analysis to the extreme statistics of the free CTRW, 
various cases of constrained CTRW can be straightforwardly studied following the analysis presented 
above. Similarly, the study of the dynamics in a disordered energy landscape generated by a CTRW 
could be done in principle from the analysis of the new fixed point of the RSRG exhibited here. 
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Appendix A. RSRG equations and fixed point for bulk and edge bonds 

Here we recall the RSRG equation for the probability distributions of barriers and lengths in a 
symmetric landscape denoted Pr(F,l) for bulk bonds and Er(F,£) for edge bonds. When convenient 
we use the notation £ = F — T. For bulk bonds it reads: 

(d r -d i )P r ((,£)= [ P r (0J 2 ) [ C dCPr(C,h)Pr(t-CJ3) , (A.l) 

Je 1 +i 2 +e 3 =i Jo 

and for edge bonds one has: 

d r E r (F,e) = -P r (0)E r (F,£) 

p pOO pOO 

+ dF i d( 3 E r (F 1 ,£ 1 )Pr(0,£ 2 )Pr(( 3> h)5(F-(F 1 +{ s )) . (A.2) 

Jii+e 2 +e3=e Jo Jo 

In terms of Laplace variable p with respect to I, one obtains: 

d r E r (F,p) = -P r (0)E r (F,p) + f dF 1 E r {F 1 ,p)P r (0,p)P r {F - F u p) , (A.3) 

Jo 

and we recall: 

Pr(C,p)=arbK CMp) (A.4) 

o(p) = sinh^/p) ' 6 r(p) = ^coth(rvp) . (A.5) 

It is straightforward to check that: 

E r (F,p) = r -i e -F^pcotHr^p) (A. 6) 

is a solution of Eq. (|A.3[) . 

Appendix B. A new fixed point for the RSRG 

We now seek a solution of (jA.ip with the scaling form: 

Pr (C,£)=T- 1 -iQr(v = f^ = f ^) , (B.1) 

which, as discussed in the text corresponds to a CTRW energy landscape with index a (we remind 
that one recovers the standard BM for a = 1). Consider Qr(ViP) the Laplace transform of Qr(r],X) 
w.r.t. A only. It satisfies the flow and fixed point equation: 

= Td r Q = Q + (1 + v)d v Q - ^-pd p Q + Q(0,p)Q(.,p) *„ Q(.,p) . (B.2) 

One can check that the solution of the fixed point equations take the form: 

Q(V,P) = a T=1 {p a )e-^=^ (B.3) 

where the functions drip) and &r(p) are given in (|A.5[) . Reexpressed in the Laplace variable p associated 
with £ we thus find: 

This is a new class of RSRG fixed points which corresponds to broad distributions of bond length: 
While there is a typical bond length, ityp ~ 

T 2 / Q the average bond length is infinite as the distribution 
does not have a first moment. Expanding at small p one has Pr(p) = 1 ~ \T 2 p a + .. from which we 
find the tail P T (£) - T 2 /£ 1+a at large £. 
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For the edge bonds one finds similarly that 

E r (F,p) = r-V^e-^^KiV*/ 2 ) (B . 6) 

is a solution (in Laplace) of Eq. (|A.3|) . Hence we see, that, as claimed in the text the CTRW fixed 
points are obtained from the BM ones (i.e. the standard RSRG fixed points) by the substitution 
p —>■ p a . 

There is however a subtlety concerning the finite size measure, i.e. the analogous of (jT]) . As 
explained in the text, it is still valid provided one interprets L as the "first arrival time" of the process 
at ul- Indeed, if one considers Zl as 

Z L = I du L f l r E T (F 1 ,£ 1 )Er(F 2 ,e 2 )S(L-(£ 1 +e 2 ))S(u L -(u + F 1 -F 2 )) 

J-oo j£ 1 ,i 2 ,F- L ,F 2 

/QQ OO n 2 — | — 1 n n 2 ~\~ 2 

du L Y, hEr{F u li) TT / p (Fj,tj) E T (F 2k+2 ,£ 2k+2 )8(L - V 4) (B.7) 

-oo k=l Jli,F x 0=2 h^Fj Je 2k + 2;F 2k+ 2 i=1 

2k+2 

x6(u L ~(u Q + ]T (-1) J+1 F,)) , (B.8) 



then one has 



(i 



where we have used lr = T 2 / a (the typical bond length) and the above forms for the fixed point 
(|B.4|) . (|B.6[) . Therefore, Zl ^ 1 for a ^ 1 and this measure is not normalized (|B.7[) if a ^ 1. The 
reason for this is that if one considers a fixed time L, one must then convolute the finite size measure 
by the waiting time function *&(£), which, in the language of the stochastic process corresponding 
to CTRW (discussed in the main text), is the probability of no jump in the time interval [£, oof. It 
is straightforward to see that its Laplace transform behaves as ^(p) = Ap a ~ 1 and the condition of 
normalization J Q e~ pL Z (p) = 1/p yields simply A = 1. 

Another way to present our result for the finite size measure is to state that the correct 
generalization of (fTJ) for models with no first moment in the bond length distribution, and statistical 
independence of successive bonds, reads: 

Z L = I du L f l r Er(F 1 ,£ 1 )Er(F 2 ,e 2 )f((£ 1 +e 2 )-L) (B.10) 

J-oo J l\M,F x ,Fi 

/QQ OO n 2 — | — 1 n n 2k~\~2 

du L J2 IrEriF^h) JJ / P{F h £j) E r (F 2k+2 ,£ 2k+2 )V(J2 £ * ~ L ) 

with the "waiting time" function ^ discussed above, such that Zl is then correctly normalized to unity. 
This measure is clearly invariant under the RSRG procedure, up to the flow of Pr and Er as described 
above. This means that the 5-function constraint on total bond length in the finite size measure is 
only possible for the BM class. It would be very interesting to study whether a starting landscape e.g. 
with a fixed total length and number of bonds will indeed flow, and in which sense, to this asymptotic 
form. At this stage our main argument is based on (i) invariance of the form (|B. 10|) under RSRG (ii) 
the CTRW interpretation given above. More generally, convergence to finite size measures has not, to 
our knowledge, been studied and is a fascinating subject. This however is beyond the scope of this 
paper. 

Appendix C. RSRG method for the extremum on the edge bond 

Let us now call Er(F, £, u, x) the joint probability that the edge bond has F, i and a maximum at x of 
value u. Then it satisfies: 

drE r (F,£,u,x) = -P r (0)E r (F, £, u, x) (C.l) 



tfrCFi ,h,u, x)P r (0, £ 2 )P r {F - F u £ 3 ) 

'max(0,r-«) ^i+^2+4i=^ 
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+ / / E r (F 1 ,t 1 ,u 1 ,x 1 )P r (0,e 2 )P r (( 3 ,£ 3 )5(F-(F 1 +( 3 )) (C.2) 

x [0(ui - (r - - ui)5(a; - Xl ) + 9(T - F 1 - Ul )5(u - (V - Fi))8(x - (h + h))] ■ 

This is equivalent to: 

drE r (F,£,u,x) = -Pr(0)E r (F,£,u,x) (C.3) 
+ / dF 1 [ 

ru px nil 

+ 6(T-u) d Ul di x I dx 1 E T {Y - uJuUuxJPr&x - £i)P r (F + u-T,£-x) . 
Jo Jo Jo 

In Laplace variables w.r.t I and x we get: 

d T E T (F, p,u,q) = -P r (0)E r (F, p,u,q)+ f dF x E T (F 1 , P ,u,q)P r (0,p)P r (F - F u p) 

Jm(0,T-u) 

+ 6{Y-u) f du 1 E r (T-u,p + q,u 1 ,q = 0)Pr{0,p + q)Pr(F + u-T,p) . 
Jo 

From the Markov property of the BM the solution of this equation must take the form: 

E r (F,e,u,x) = A r (F,u,x)B r (u + F,£-x) , (C.4) 

E r (F,p,u,q) = A r (F,u,p + q)B r (u + F,p) , (C.5) 

where A is the sum over paths starting at 0, ending at x, u constrained to remain on interval [— F, u] 
and with no return of more than T, while B are paths starting at x, u ending at £, —F constrained to 
remain on the interval [-F, u] with no return of more than T. For T > u the return constraint does 
not play a role for A, for T > u + F it does not play a role for B. Whenever return constraint does 
not play a role one has: 

MF,u, P )= .^f^ , (C.6) 
smh(y/p(u + F)j 

B r (u + F, P )= ^ (C.7) 
smh{y/p(u + F)) 

i.e. independent of T. For smaller T one must solve the full equation above with initial condition 

Er(F = T,£,u,x) = S(u)S(x)Er(F = T,£). We will not attempt this here as we are mostly interested 
in the large T limit. 

Hence our final result for the solution at large T reads: 

W fF ii ii g\ — sinhj^p + qF) VP (r „, 

*r[*,P,u,q) Ts - mHv¥T ^ {u + F))sinHV p {u + F) y W 

One can check that 

1 
duE r (F,p,u,0) = -e~^ F , (C.9) 
o 1 

which is the correct result. 

We can now compute E r (F,£,u,x) by performing a double inverse Laplace transform. One has 

indeed 

x ( sinh(ypT^) \ = A M^n) ^ x 

L1 ^Unh(Vm(u + F))J ~ ^o ] (u + F) 2 ■ (U1Uj 

Note that this yields the identity (setting p = and taking the Laplace transform wrt to x of both 
sides) 

2,^^ I - sinh ^ F > (C 11) 



L 
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which we also checked with Mathematica. 
Next we have 

Note that this yields trivially the identity (setting x — and taking the Laplace transform wrt to y of 
both sides) 

2tt \ 1 ^ V (-l) m+1 = ^ - , (C.13) 

(u + F)3^ J p+^f^ sinh(Vp(u + F)) 1 ' 

which we also checked with Mathematica. Note that the sum in the left hand side has to be understood 

as 

oo 2 i 00 2 

27T 2 , „ V (-1)" 1+1 V ^^= lim 2^ 2 * y a m+1 . (C.14) 

l U + ^ J m=0 P + («+-F) 2 1 ' ' m=0 P + (u+FY 

Finally, combining the both Laplace inversion in Eq. (IC 121 IC.14[) one obtains the formula (jTTJ) 
given in the text. 

Appendix D. Details about reflected BM, excursions and meanders 

Here we give the details about the computation of the joint distribution Pi,(u m , x m , u, x, ul) of the 
global minimum u and maximum u m and their positions x and x m . We start with the expression given 
in the text in Eq. |20|) : 

P L (u m ,x m ,u,x,u L ) = 6[x m - x) I du 2 I dx 2 P L (u,x,u m ,x m ,u 2 ,x 2 ,u L ) (D.l) 



+ 6{x — x m ) \ du\ I dXlPL(ui,Xl,Um,%m, u,x,u£) , (D.2) 
Ju JO 

Let us first focus on the first term. The integral over the space variable x 2 can be done by noticing 
that it can be written, formally, as 

/ du 2 I dx 2 P L (u,x,u m ,x m ,u 2 ,x 2 ,u L ) = / du 2 J dx 2 F(L - x 2 )G(x 2 - x m ) 

Ju J £ m J u J x m 

du 2 / dy 2 F(L - x m - y 2 )G(y 2 ) , (D.3) 



where the functions F, G can be read straightforwardly on Eq. Q19p . The integral over y 2 can be 
performed by taking the Laplace transform with respect to L = L — x m . This yields 

dye~ p£ [ " du 2 [ dy 2 F{L - y 2 )G{y 2 ) (D.4) 



Ju JO 

oo 



V 47rVm 2 (-l)^+ m i Sm( "-^^ (D.5) 



X 



ul 



OO rj.y. / 7T (u rn U l ) 



du 2 £ (47r 3 )( _ 1) » 2+m2n2m 2 Sln u m -u 2 "2; 1 

«2,m 2 =0 V ' P+ (u m -u 2 ) 2 P+ (u m -u 2 y> 

Now one can use the identities in Eq. (jC.lli IC.13[) to perform the sums over n\ , m,\ to obtain 

dye- pi { " du 2 [ dy 2 F{L - y 2 )G(y 2 ) (D.7) 



Ju JO 
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= Y 47r 3 n 1 m?(-l)" 1 + m 1 ^S^e-^^ ( " ?a;+m?(a;m - a;)) 
x f L du2 VP^f^-u L )) (D8) 



u sinh (y/p(u m - u 2 )) 



sinh (^/p(u L - u)) 



= Y 47r 3 n 1 m?(-l)^^^^^e-^^ (n ? a+mi(xm - x)) - 
One can then inverse the Laplace transform (see Eq. (|C.12j) ) to obtain 



(D.9) 



du 2 / dX2P L (u,X,U m ,%m,U2,X2,U L ) (D.10) 

y 47r 3 n 1 m^(-l)" 1 + mi Sm( "" l) e '(^^ ( "' 3:+m ' (3: '"'" )) (D.ll) 



y 2 7 r(-l)-+ 1 n 2 ^S^e-<^ (i -- ) . (D.12) 



' n(u L —u) \ 22 
\— ....... ' 

X 

n2— 

The second term in Eq. ([20]) can be computed in a similar way to get finally the expression given 
in Eq. (j2"2"j) . 

Appendix E. Marginal distribution of x m for the Bessel bridge : link with the result of 
Pitman and Yor 

We start from the joint distribution P^ cssclblldg0 (u m , x m ) given in the text in Eq. ([50]) . We remind it 
here: 

pBessel bridge/ \ _ olj \ " J _^m ^— „a /-p,,-* 

r(l + v)uU 2u Jv+l(jv,n)J»+l(jv,m) 

The marginal distribution of the maximum u m is obtained after integration over x m . It reads 

pBessel bridge/ \ _J_p Bessel bridge ( u m \ 

■2 -2 
, 1 ,1 K m 

o AV+lAV+l _ o 2^ 

pBessel bridge/ \ _ 2 v ' n Jv ^ m x " (T? o\ 

1 ' T(l + ^+2- 2-J J v+1 (j V)n ) J v+1 0V,m) # m - Jin ' ^ ■ ' 

while Pitman and Yor found in Ref. [3U] (note that the value of the diffusion constant used in their 
paper is smaller than ours by a factor of 2) 

P ^^ {x) = ^(^_^^ )eK J^\\ . (E . 3) 



We want to show that these two expressions in Eq. (IE.2P and Eq. (|E.3|) are identical and therefore we 
want to show the identity 

■ 2 -2 

V—' Jv.n Jv.m e x e an 

r 3+2i/ 7 ~ (j 1 / ~ (j \ j2 ~2 V ■ ' 

9 ,-2i/+2 / -2 \ 9 , 9 oo . 2 j/ / -2 \ 

z Ji/,n / JiAn \ T ^ v Jiy.n / Ji^.n I /-p _x 

" ^T5 2_, J2 +I (^ n ) J ~ ~^J+3 2_, j2 +i{jun) «P J ■ (E.5) 

Now in the double sum over n, m in the left hand side of Eq. (|E.4j) we separate out the m = n term 
since it exactly cancels the first term in the right hand side of Eq. (|E.4[) . Finally using a formula 
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obtained by Pitman and Yor in Ref. [3D] (see their formula (125), notice however that there is a 
misprint in their formula (124)) one obtains 

-2 -2 

2 L j^( 7 U^f7 ) 7 2 - 7 2 --2(^+1)^ exp --^ , (E.6) 

which shows finally the identity in Eq. (1E.4|) . 
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